Evidence for Dynamic Charge Stripes in the Phonons of Optimally Doped YBCO 
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Inelastic neutron scattering investigations on optimally doped YBCO revealed a very pronounced 
temperature dependence in the Cu-O in-plane bond-stretching vibrations along the (010) direction: 
a shift of spectral weight by at least 10 meV has been observed in a narrow range of wave vectors 
halfway to the zone boundary. The temperature evolution starts at about 200 K, well above the 
superconducting transition temperature. The displacement pattern of the anomalous phonons indi- 
cates a dynamic one-dimensional charge-density modulation within the planes, but excludes a direct 
contribution from the Cu-O chains. 

PACS numbers: 74.25.Kc, 63.20.Bk, 74.25.-Jb, 74.72. Bk 



There is widespread belief that the electronic struc- 
ture of high temperature superconductors is character- 
ized by an inhomogeneous distribution of charge and 
spin in which the spins on the copper atoms form an- 
tiferromagnetic stripes, separated by domain walls car- 
rying the charge. Although theory P, Q suggests that 
stripe formation should be a universal property of the 
cuprates, solid experimental evidence for stripes has 
been found so far only for a particular compound, i.e. 
La1.45Ndo.4Sro 15CUO4 3], in which the stripes are static 
and apparently detrimental to superconductivity. In su- 
perconducting compounds, the stripes are expected to be 
dynamic in nature, but dynamic stripes are much more 
difhcult to detect and so the experimental basis in favor 
of dynamic stripes is far from being solid. Inelastic neu- 
tron scattering can be used to observe the dynamic spin 
structure by the interaction between the neutron spin and 
the electron spin. Magnetic fluctuations consistent with 
the stripe phase picture have indeed been observed in 
the system La2-3;Sr3;Cu04 0,0 whereas the interpreta- 
tion of the results for YBa2Cu307_a: remains controver- 
sial Inelastic neutron scattering can be used to look 
for fingerprints of dynamic charge stripes as well, but in 
this case, evidence for stripes is obtained only indirectly 
by searching for effects of the dynamic charge structure 
on the atomic vibrations. That is to say, phonons with 
a displacement pattern closely related to that resulting 
from the inhomogeneous charge distribution are expected 
to show an anomalous behavior somewhat reminiscent of 
the phonon anomalies observed in one-dimensional con- 
ductors at temperatures above the Peierls transition, i.e. 
as a precursor phenomenon to charge-density-wave for- 
mation p^. 

In this Letter, we report observations by inelastic 
neutron scattering on optimally doped YBa2Cu3706.95 
which fit nicely to the predictions of recent theoretical 
work 0, 13 as to phonon anomalies resulting from an 
interaction with low- lying collective charge excitations. 



We show that a pronounced downward shift of spectral 
weight occurs for Cu-O in-plane bond-stretching vibra- 
tions propagating along the b-direction within a narrow 
range of wave vectors on cooling below room tempera- 
ture. This result appears to be difficult to explain other 
than with the development of dynamic charge stripe 
correlations. We note that first results on an unusual 
temperature dependence of phonons in optimally doped 
YBCO were already reported in 0. An important step 
for achieving a better understanding of these effects [Tj 
has been to study the phonon branches of A4-symmetry 
instead of those of A 1 -symmetry : these branches dif- 
fer from each other in that the elongations in the Cu-O 
bi-layer are in-phase for Ai-symmetry and out-of-phase 
for A4-symmetry. We found that the behavior of the 
two types of bond-stretching branches is in principle very 
similar. However, anticrossings of the in-plane-polarized 
bond-stretching branch with c-axis- polarized branches 
leads to a formidable complication in the case of the 
branches of Ai-symmetry 



npiicatic 



Inelastic neutron scattering measurements were per- 
formed on a composite sample of optimally doped 
YBa2Cu3 6.95 (Tc — 93 K) consisting of three twinned 
single crystals with a total volume of 1.5 cm'^ . De- 
tails of the sample are given elsewhere The neu- 
tron experiments were carried out on the IT triple-axis 
spectrometer at the ORPHEE reactor of the Laboratoire 
Leon Brillouin at Saclay, France. The (220) reflection 
of a Cu crystal was used to monochromatize the inci- 
dent neutrons in order to achieve high resolution. Py- 
rolytic graphite (PC) (002) was used as analyzer crystal. 
Both crystals were horizontally and vertically focusing. 
A PC fllter was placed into the final beam to suppress 
higher order contaminations. The measurements were 
performed at various temperatures between T=10 K and 
300 K. The experiments focused on the Cu-O in-plane 
bond-stretching vibrations in the (100) and in the (010) 
directions because it was known from previous investiga- 
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FIG. 1: Dispersion relation of the out-of-phase Cu-0 bond- 
stretching vibrations in optimally doped YBCO as seen on a 
twinned sample. The energies are plotted against the nominal 
wave vector, which is calculated from an average lattice con- 
stant aai, = (a+b)/2. The actual wavevectors for the scattering 
from each of the two twin domains differ by 2(5=27r/a-27r/b. 
The arrows show the zone center for phonons propagating 
along the a or the b direction, respectively. The ellipsoid 
depicts the instrumental resolution. 



tions [1^ [13, Hal that these branches show an anomalous 
dispersion at low temperatures. 

The twinning of our sample was, of course, a diffi- 
culty for the assignment of the neutron peaks to the 
a or the b direction. This problem was overcome by 
the technique described in refs. 9,15, i.e. by exploit- 
ing the fact the dispersion relation for phonons propa- 
gating along (100) is shifted in reciprocal space against 
that of phonons propagating along (010) due to the dif- 
ference in lattice constants of nearly 2% (see Fig. 1). 
A further strong support for our assignment came from 
previous measurements on a detwinned sample |l2l Il3j | . 
Unfortunately, this sample cleaved after the first cool- 
ing cycle and could therefore not be used for the tem- 
perature dependent measurements. As can be inferred 
from Fig. 1, measurements made for momentum trans- 
fers Q=(3,0,2)-|-q or Q=(3,0,2)-q, respectively, have both 
their pros and cons: for Q=(3,0,2)-q, focusing leads to 
relatively narrow phonon lines, but the peaks for the a* 
and the b* directions are rather close in energy, whereas 
the opposite is true for Q=(3,0,2)-|-q. Therefore, both 
configurations have been tried to elucidate the behavior 
of the bond-stretching phonons. 

A major result of our study is that strong temperature 
effects were observed only for b-axis polarized phonons 
and only in the vicinity of the wave vector q=(0, 0.27,0). 
This led us to carry out a detailed temperature study 
for this particular Q. We found a pronounced shift of 
spectral weight from energies around 59 meV to energies 
around 49 meV when lowering the temperature from 300 
K to 10 K (Fig. 2a,b). The intensity loss at E«59 meV 
appears to be not fully compensated by the gain at ERi49 
meV. Therefore, we made further scans covering a larger 
energy interval to search for the missing intensity. The 
data suggest that some spectral weight is shifted down 
to E«37 meV (Fig. 2c) with the caveat that it is dif- 
ficult to correct for the temperature dependence of the 



background for such wide energy intervals. 

Fig. 2b shows that the temperature evolution starts at 
temperatures far above Tc and also continues on further 
cooling below Tc. The question as to whether there is 
a subtle influence of Tc on the temperature dependence 
of the phonon under discussion will have to be addressed 
in a future investigation using finer steps in temperature. 
Fig. 2b further suggests that the phonon softening pro- 
ceeds by a shift of spectral weight and not a peak shift, 
i.e. a gradual depletion of the phonon intensity observed 
around 59 meV at T=300 K and a simultaneous build- 
up of intensity at much lower energies. This issue was 
further explored by scans at the nominal wave vector 
Q=(2.7,0,2), corresponding to the same q along b* as 
at Q=(3.25,0,2) because focusing leads here to a better 
energy resolution. Unfortunately, the contributions from 
the a* and the b* directions show up at nearly the same 
energy at T=200 K. On cooling, a pronounced loss in 
intensity is observed which corresponds very well to that 
observed at Q=(3.25,0,2). However, no peak broadening 
was seen at any temperature demonstrating that all the 
lost intensity is shifted out of the energy window of these 
scans except for a modest intensity gain at the lowest end 
of the window (50 meV) (Fig. 3b). We note that this in- 
tensity loss rapidly disappears when going to smaller q 
(Fig. 3a) . On the high q side, the temperature effect be- 
comes rapidly weaker as well and leads to a softening of 
the bond-stretching phonon peak of a few meV only. At 
the zone boundary, there was no discernible temperature 
dependence of the bond-stretching mode (Fig. 3c). 

The phonon anomaly described above bears certainly 
some resemblance to the Kohn anomalies observed in 
one-dimensional conductors |17| . i.e. precursor phenom- 
ena to a charge- density- wave transition occurring at low 
temperatures. Therefore, one might conjecture that it 
is related to flat parts of the Fermi surface related to 
the oxygen chains However, we emphasize that the 
phonon anomaly is observed in a branch of A4-symmetry 
for which the chain O elongations are zero. This fact ex- 
cludes a trivial relationship of the anomaly with chain 
O related electronic states. Having ruled out a direct 
association with the Cu-0 chains, we conclude that the 
observed anomaly is evidence for a significant electronic 
anisotropy, consistent with nematic behavior p^ . 

Could the phonon anomaly be the result of Fermi sur- 
face nesting within the Cu-0 planes ? A check of this 
possibility is provided by recent results of density func- 
tional theory [13 . These calculations indeed predict a lo- 
cal frequency minimum in a branch with bond-stretching 
character at qRi(0, 0.25,0); however, this branch is of 
A 1 -symmetry and involves primarily chain O elonga- 
tions. This feature is not unexpected in view of the one- 
dimcnsional nature of the Cu-0 chains and may well be 
related to static charge modulations detected within the 
Cu-0 chains by scanning tunneling microscopy (STM) 
[25I |2^ . On the other hand, the calculations do not 
predict a similar anomaly in a branch of A4-symmetry. 
Rather the theoretical results are in surprisingly good 
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FIG. 2: (a) Energy scan taken at a nominal momentum trans- 
fer of Q=(3.25,0,2). Due to the twinning of the sample this 
Q corresponds to a reduced wave vector q=0.22 along a* or 
q=0.28 along b*, respectively. The 300 K data have been 
corrected for the increase in multi-phonon scattering on rais- 
ing the temperature from T=10 K. (b) Smoothed intensity 
difference with respect to the intensity measured at T=300 
K. (c) Difference between the intensities measured at T=10 
K and T=200 K, respectively, as measured in an independent 
run over a larger energy interval. 



agreement with the 200 K data (Fig. 4). The authors 
of that study were concerned that the high effective elec- 
tron temperature (0.2 eV) used in their calculations to 
improve convergence might have suppressed the anoma- 
lous softening at q=0.25. Therefore, additional calcula- 
tions were made using a much smaller effective electron 
temperature (0.02 eV). However, this had barely any in- 
fluence on the calculated phonon frequencies [23|. Ap- 
parently, the anomalous low-temperature dispersion of 
the in-plane bond-stretching phonons cannot be under- 
stood within the framework of such a theory. For this 
reason, it seems unlikely that the observed anomaly can 
be attributed to a Fermi surface nesting effect. 

A better candidate to explain the observed anomalies 
is the coupling of the phonons to charge stripes, similar to 
what has been proposed in recent theoretical work 0| . 
As indicated in Fig. 5, the displacement pattern of the 
anomalous phonons is such as to favor dynamic charge 
accumulation on (approximately) every fourth row of 
atoms along b. The idea of charge stripe formation moti- 
vated us to search for elastic superlattice peaks related to 
static charge density wave formation, but we found none 
|25| . This means that if our interpretation of the phonon 
anomaly is correct, the charge stripes in optimally doped 
YBCO are really dynamic in nature. 

The above reasoning assumes that the charge modula- 
tion is along b, which means that the stripes must run 
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FIG. 3: Energy scans taken at different temperatures. Lines 
correspond to fits with one or two Gaussians. The 100 K data 
points have been omitted for clarity. Due to the twinning of 
the sample, the nominal wavevectors indicated in the figure 
correspond to q=0.23 a* and q=0.17 b* (a) or q=0.33 a* 
and q=0.27 b* (b) or q=0.47 a* and q=0.53 (=0.47) b* (c), 
respectively. The double peak in (a) corresponds to the a axis 
mode (67 meV) and the b axis mode (63 meV), respectively. 
These two peaks merge to a single one in (b) but are again 
well separated in energy in (c). The 54 meV peak in (c) 
also contains contributions of the bond-bending mode (see 
the dispersion relation plotted in Fig. 4). 



parallel to a, an orientation that would be orthogonal to 
that implied by studies of magnetic fluctuations in un- 
derdoped YBCO [Ull^l. On the other hand, a charge- 
density-wave instability (CDW) along charge stripes has 
been discussed theoretically in j2J|| : formation of metallic 
stripes would make the electronic properties of a cuprate 
similar to that of one-dimensional conductors and hence 
could lead to the same type of transition which is typ- 
ically observed in such materials. Of course, this ar- 
gument tacitly assumes that the onset temperature for 
stripe correlations is quite high which is an unresolved 
issue. 

Comparing the data for optimally doped YBCO (07) 
with those of underdoped 06.6 we find that the low 
T dispersion of 06.6 definitely shows an anomaly similar 
to that observed in 07 although not quite as clearly. Re- 
inspection of temperature-dependent data also revealed 
a significant temperature effect which was overlooked in 
[27I I because it does not show up as a peak shift but again 
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FIG. 4: Comparison of the phonon dispersion of the Cu-0 in- 
plane bond-stretching phonons (E>55 meV) and of the Cu- 
O bond-bending phonons (E<55 meV) of A4-symmetry as 
calculated from density functional theory by [isll with the 
experimental data obtained at T=200 K. The red line is a 
guide to the eye. 
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as a shift of spectral weight. This indicates that the 
phonon anomaly reported in this paper is not restricted 
to optimally doped samples. 



In summary, we observed a pronounced temperature 
effect in the Cu-0 in-plane bond-stretching vibrations 
strongly indicative of dynamic charge stripe formation 
in optimally doped YBCO. This result certainly adds 
credence to long-held theoretical expectations that the 
cuprates should form a charge-inhomogeneous state. On 
the other hand, several questions remain to be answered: 
Why is the effect, if present at all, much less pronounced 
along the a direction ? Why is the effect not seen in op- 
timally doped LSCO ? Last but not least, how does 
one reconcile the magnetic response of optimally doped 
YBCO with the idea of a dynamic stripe phase ? As will 
be discussed in a separate paper [2^ antiferromagnetic 
spin fluctuations have indeed been observed on the same 
sample, but their dispersive nature and their apparent 
in-plane isotropy do not fit to the naive expectations of 
the classical stripe picture ■ 



FIG. 5: Displacement pattern of the longitudinal Cu-0 in- 
plane bond-stretching phonon with wave vector q=(0, 0.25,0). 
Only the elongations within one Cu-O plane are shown. The 
elongations in the other layer of the bi-layer are in-phase for 
phonons of Ai-symmetry and out-of-phase for phonons of A4- 
symmetry, respectively. Circles and full points denote Cu 
atoms and O atoms, respectively. Filling of circles indicates 
dynamic charge accumulation on every fourth row of atoms. 
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